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CLT Abstract 

^ ! 

We calculate the Next-to-Leading Order (NLO) QCD corrections to Zbb production in hadronic 

. collisions including full bottom-quark mass effects. We present results for the total cross section 

>. 

and the invariant mass distribution of the bottom-quark jet pair at the Fermilab Tevatron pp 

o' 

collider. We perform a detailed comparison with a calculation that considers massless bottom 
\Q ■ quarks, as implemented in the Monte Carlo program MCFM. We find that neglecting bottom- 

o, 

quark mass effects overestimates the total NLO QCD cross section for Zbb production at the 

o 

Tevatron by about 7%, independent of the choice of the renormalization and factorization scales. 
^ Moreover, bottom-quark mass effects can impact the shape of the bottom-quark pair invariant 

mass distribution, in particular in the low invariant mass region. 
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I. INTRODUCTION 



One of the main goals of high-energy collider experiments is the elucidation of the mecha- 
nism of Electroweak Symmetry Breaking (EWSB) as well as the exploration of energy scales 
beyond the weak scale, where physics beyond the Standard Model (SM) is expected. The 
hadronic production of weak gauge bosons in association with a bottom-quark pair plays 
a crucial role in some of the current studies of EWSB and beyond the SM physics at the 
Fermilab Tevatron pp collider [l[ [2], . Wbb and Zbb production processes represent 

the major irreducible backgrounds to the main search modes for a light SM-like Higgs boson 
at the Tevatron, i.e. WH and ZH with H — > bb. Wbb also accounts for one of the most 
important backgrounds to single-top production, pp tb, tb with t(t) — > Wb(b), which tests 



the fundamental structure of the Wtb vertex at the Tevatron 

say , llOj. Finally, Zbb is 

a background to searches for Higgs bosons in models with enhanced bottom-quark Yukawa 
couplings, such as the Minimal Supersymmetric Standard Model (MSSM) with large tan j3, 
where Hbb with H — > yU + /i~,r + r~ is an interesting discovery channel 11]. 

The hadronic cross sections for W/ZH associated production have been calculated includ- 
ing up to Next-to-Next-to-Leading Order (NNLO) QCD corrections 3, [l3, Q and 0(a) 



electroweak corrections [15j. Single-top production has been calculated at Next-to-Leading 
(NLO) in QCD [l6, 17, 18, 19, 2fl 21, 22, 23, 



ft 

123 



241 ]. and including one-loop electroweak (SM 



and MSSM) corrections [25|], while the cross section for Hbb associatedproduction is known 



including NLO QCD corrections and full bottom-quark mass effects 26 



27 



28 



29]. 



To fully exploit the Tevatron's potential to detect the SM Higgs boson or to impose 
limits on its mass, it is crucial that the dominant background processes are also under 
good theoretical control. In the present experimental analyses 1 , the effects of NLO QCD 
corrections on the total cross section and the dijet invariant mass distribution of the W/Z bb 
background processes have been taken into account by using the Monte Carlo program 
MCFM |ao|. In MCFM, the NLO QCD predictions of both total and differential cross 
sections for the se p rocesses have been calculated in the zero bottom-quark mass (m^ = 0) 



approximation 



31 



32 



33| . using the analytical results of Refs. [34j, [35(. From a study of the 



Leading Order (LO) cross section, finite bottom-quark mass effects are expected to affect 



1 For updated results, see the CDF and £>0 websites at http://www-cdf.fnal.gov/physics/exotic/exotic.html 



and |http: / / www-dO.fnal.gov /Run2Physics /WWW/results /higgs. htm 
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both the total and differential W/Zbb cross sections, mostly in the region of small 66-pair 



invariant masses 



33j. Indeed, since this kinematic region of small 66-pair invariant masses 



contributes considerably to W/Z + nj production (n = 1, 2), where at least one o 



the jets 



is a 6-jet, bottom-quark mass effects cannot be neglected as discussed in Refs. |36l . [37| (for 
n = 2) and in Ref. [38| (for n = 1). Given the variety of experimental analyses involved in 
the search for W/ZH associated production, single-top and Hbb production, it is important 
to assess precisely the impact of a finite bottom-quark mass over the entire kinematical reach 
of the process, including the complete NLO QCD corrections. 



In Ref. 



39( we have performed a study of NLO QCD cross sections and invariant mass 



distributions of the bottom quark pair in Wbb production at the Tevatron, including full 
bottom-quark mass effects. We found that bottom-quark mass effects amount to about 8% 
of the total NLO cross section at the Tevatron and are mostly visible in the region of low 
66-pair invariant mass. 

In this paper, we compute the NLO QCD corrections to Zbb hadronic production, includ- 
ing the full bottom-quark mass effects. We consider all partonic processes that contribute 
at 0(aal), i.e. NLO QCD corrections to qq — > Zbb and gg — *> Zbb and the tree-level process 
q{q)9 ~^ Zbbq(q). We present numerical results for the total cross section and the invariant 
mass distribution of the bb pair for the Tevatron pp collider, including kinematic cuts and a 
jet-finding algorithm. In particular, we apply the kr jet algorithm and require two tagged 
6-jets in the final state. Using the MCFM package [30], we compare our results with the 
corresponding results obtained in the mj = limit. Numerical results for both Zbb and Wbb 
production at the Large Hadron Collider will be presented in a separate publication [iol ]. 

The outline of the paper is as follows: in Section [TT] we briefly discuss the technical details 
of our calculation, while we present numerical results and a discussion of the bottom-quark 
mass effects in Section IIHI Section [IV] contains our conclusions. 
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FIG. 1: Tree-level Feynman diagrams for qq — ► Zbb. The circled crosses indicate all possible 
insertions of the final-state Z boson leg, each insertion corresponding to a different diagram. 



FIG. 2: Tree-level Feynman diagrams for gg — ► Zbb. The circled crosses indicate all possible 
insertions of the final-state Z boson leg, each insertion corresponding to a different diagram. 

II. CALCULATION 
A. Basics 

The total cross section for pp(pp) — > Zbb at O(ao z s ) can be written as follows: 



where denote the parton distribution functions (PDFs) for parton % in a proton (an- 

tiproton), defined at a generic factorization scale fj,f — fx. The sum runs over all relevant 
subprocesses contributing to the hadronic cross section initiated by partons i and j. The 
partonic cross section for the subprocess ij — > Zbb (+k) is denoted by <J^ LO and is renor- 
malized at an arbitrary scale /i r = //. If not specified otherwise, we assume the factorization 
and renormalization scales to be equal, fif = fi r = fi. The factor in front of the integral is a 
symmetry factor that accounts for the presence of identical particles in the initial state of a 
given subprocess (Sij is the Kronecker delta). The partonic center-of-mass energy squared, 





[jVfarijfr^rio^iz^ri + (a* - x 2 )] , (1) 
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s, is given in terms of the hadronic center of mass energy squared, sh, by s = x±x 2 sh- 
The NLO QCD partonic cross section reads: 



a 



NLO, 



(x 1 ,x 2 ,/J l ) = &lf > (x 1 ,x 2 ,n) + 5crf j LO (x 1 ,x 2 ,fi) 



NLO/ 
ij 



(2) 



where a\j (xx,x 2 ,ii) denotes the 0(aa 2 s ) LO partonic cross section and (x±, x 2 , p) 

describes the 0(a s ) corrections to aff(xi, x 2 , p). The LO pp(pp) — > Zbb process receives 
contributions from qq and gg initiated processes, as shown in Fig. [1] and Fig. [2j respectively. 
The NLO QCD corrections, 5afj LO , receive contributions from qq, gg, qg and qg initiated 
processes and can be decomposed in the following way: 



5a 



NLO 

ij 



j d {PS-,) J^l Airt(ij - ^)| 2 + J d (PS A ) ^|Acai(u -> Zbb + fc)| 5 



= cr. 



virt 



a 



real 



(3) 



where the term integrated over the phase space measure d (PS 3) corresponds to the virtual 
one-loop corrections with three particles in the final state, while the one integrated over 
the phase space measure d (PS 4) corresponds to the real tree-level corrections with one 
additional emitted parton. The sum ^ indicates that the corresponding amplitudes squared, 
|^virt(reai)(u ~~ > Zbb(+k))\ 2 , have been averaged over the initial-state degrees of freedom 
and summed over the final-state ones. The phase space integration has been performed 
using Monte Carlo (MC) techniques using the adaptive mult i- dimensional integration routine 



VEGAS 



4l|. 



We have improved on the massless calculation of Refs. [32|, l33( by considering a fully 
massive 6-quark both at the level of the scattering amplitude and in the integration over the 
phase space of the final-state particles. We keep the Z boson on-shell, though the extension 
to include its leptonic decays does not present in principle any special complications. Because 
of the complexity of this calculation, all results have been cross-checked with at least two 
independent sets of codes. The analytic calculation of the scattering amplitudes has been 
implemented using, at different stages, FORM j^J, TRACER 43], Mathematica and Maple. 



Final numerical resu 



ts have been obtained with codes built in C and FORTRAN, and we 



have used the FF [44| and Madgraph [45|, |46|, |47j packages for cross-checks. 

The 0(a s ) corrections to Zbb production are similar in structure to the NLO QCD 
corrections to Wbb production [39] and to Higgs production in association with top quarks 



(Hit) 48|, |49j. We therefore only summarize below the most important features of the 



calculation and refer to the Wbb and Htt papers cited above and to Ref. 50j for more 
details. The 0(a s ) corrections to qq — > Zbb can be derived from the 0(a s ) corrections 



48] (with t <r+ b and 



to qq' — ■> Wbb production [39J (with W <-> Z) and to gg — >■ ifi? 
H Z), while the C?(a s ) corrections to the initiated production process can be 
obtained from the 0(a s ) corrections to gg — > iftt [49] (with £ «-> 6 and H Z). The 
gg, gg initiated processes appear in both and iftt NLO QCD calculations with W <-> Z 
and H ^ Z,t ^ b, respectively. Note that, when applying the results of and iiftt 
production of Ref. and Refs. Q, 49], one also needs to replace respectively the Wff 
and fermion Yukawa couplings by the V — A coupling of fermions to the Z boson: 




2 sin 9 W cos 9 W 



(4) 



where 6 W is the weak mixing angle and the vector, g } Vl and axial- vector, g } Al couplings for 
the Zff vertex are given by: 

g f v = Ti -2 sin 2 9 w Q f , g{ = , (5) 

with Tl denoting the third component of the weak isospin and Q / the electric charge of the 
fermion /. Moreover, as we are considering an on-shell Z boson, we have summed over its 
polarizations as follows: 

^<?M<r(pz) = -ir + (6) 

where M z is the mass of the Z boson. 



B. The virtual cross section <5-y„! rt 



The 0(a s ) virtual corrections to the partonic tree-level gg — > Zbb and gg — > Zbb produc- 
tion processes consist of self-energy, yertex, box and pentagon diagrams, as shown, for the 



Htt-like part, in Figs. 2-4 o 



Ref. 



481 ] and Figs. 2-5 of Ref. [49j], respectively (for a full set 



of diagrams see also Ref. |50|). The contributions to aj" 1 in Eq. <^ can be written as: 



EE 

D 



D 



(7) 



where Aq is the tree level amplitude and Ad denotes the amplitude for the one-loop diagram 
D, with D running over all self-energy, vertex, box and pentagon diagrams corresponding 
to the r/'-initiated subprocess. 

The calculation of each virtual diagram (Ad) is performed in the same way as outlined 



in Refs. 



48l . |49| . i.e. Ad is calculated as a linear combination of Dirac structures with 



coefficients that depend on both scalar and tensor one-loop Feynman integrals with up to 
five denominators. We solve the one-loop integrals in the coefficients either at the level 
of the amplitude or at the level of the amplitude squared (see Eq. (J7|)). These two in- 
dependent approaches allow us to thoroughly cross-check the calculation of each individual 
diagram. Indeed, the tensor structures present in the one-loop integrals of the amplitude are 
typically different from the ones present in the amplitude squared, as one can perform non- 
trivial reductions of the latter by canceling dot-products of the integration momentum in 
the numerator with denominators in the Feynman integrals. In this way, the final analytical 
expression of a given diagram ends up being represented in terms of different building blocks. 
A possible incorrect relation between the building blocks would then naturally produce a 
discrepancy between the two approaches. 

Tensor and scalar one-loop integrals are treated as follows. Using the Passarino-Veltman 
(PV) method 5l|, |52[, the tensor integrals are expressed as a linear combination of tensor 
structures and coefficients, where the tensor structures depend on the external momenta and 
the metric tensor, while the coefficients depend on scalar integrals, kinematics invariants 
and the dimension of the integral. Numerical stability issues may arise at this level as a 
consequence of the proportionality of the tensor integral coefficients to powers of inverse 
Gram Determinants (GDs), as discussed in detail in Ref. [49|] , although for Zbb production 
the problem is considerably more serious. 

These numerical instabilities can be considered as "spurious" or "unphysical" divergences, 
since it is well known that only two-particle invariants can give rise to a physical singularity. 
Indeed, these spurious divergences cancel when large sets of diagrams are combined [35(, such 
as, for example, when one combines gauge invariant sets of color amplitudes (i.e. amplitudes 
with a common color factor). As we have expressed our calculation in terms of invariants, 
and we employ a standard basis of scalar integrals, the full cancellation only occurs between 
numerator and denominator at the numerical level, often between fairly large expressions. 

For this reason we have chosen to organize the diagrams, at certain stages, into gauge in- 
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variant color amplitudes, that is, into coefficients of the same color structure (see Ref. [50| ) . 
This allows a better handling of the spurious singularities and a natural way to make in- 
ternal cross-checks and cross-checks with new techniques. When we consider these gauge 
invariant sets of color amplitudes and full analytical reductions of all tensor integrals, we find 
cancellation of some powers of GDs, which improves the numerical stability of our code, so 
that when integrating over the Zbb phase space, using MC techniques, we obtain statistical 
errors below 0.1% for total cross sections. 

The fully reduced numerical codes are often more demanding computationally, and be- 
cause of that we have built master codes that use them only when close to regions of phase 
space where certain problematic GDs become small. All this is found particularly useful 
when considering higher rank D-PV functions (we have up to DA-PV functions in our cal- 
culation) as well as E-PV functions. This technique would probably break down if one were 
to extend it to processes with even more legs, and the use of helicity amplitudes would in 
this case be preferable. 

In the case of pentagon diagrams, a powerful and convenient check consists of reducing 
consistently all E-PV functions by canceling systematically, at the level of the amplitude 
squared in Eq. (j7]), all possible vector products containing the loop momentum in the numer- 
ator with some denominators. This is possible as, in the pentagon topology of our process, 
each leg has an outgoing momentum which is on-shell, corresponding basically to one of 
the external initial or final particles of the subprocess. One then ends with expressions for 
each pentagon diagram containing purely scalar pentagon integrals, or tensor integrals with 
fewer than five denominators, improving considerably the numerical stability. We compared 
analytically these reductions to the non-reduced expressions by using the full reduction of 
all tensor integrals to scalar integrals, and found agreement. 

We also checked parts o 



our result by using unitarity techniques [35|, specifically the 
quadruple-cut technique |53j]. As shown by Britto, Cachazo and Feng (BCF), from any set 
of Feynman diagrams (or more generally from any tensor integral 54J) one can extract the 
coefficient of a given scalar box integral by cutting the four corresponding propagators (see 
Fig. [3]), i.e. by replacing i/(p 2 — m 2 + ie) — > 2-k5^ + \p 2 — m 2 ) for each cutted propagator of 
momentum p and mass m. This effectively freezes the momentum integration, and replaces 
it by a set of algebraic equations which determine the loop momentum entirely. We solved 
this set of equations by using a BCF ansatz [53] , and then compared the result to the corre- 



S 




FIG. 3: Quadruple cut 



53| check of the calculation of a box diagram involving a top-quark loop. 



It corresponds to two Feynman diagrams given by the two possible orientations of the fermion line. 



sponding box coefficient extracted from our analytic expression, and found agreement (for 
more details and specific solutions for the topology in Fig. [3] see Ref. 50|). This is a rather 
non-trivial check for the set of E-PV and D-PV functions we have employed at different 
stages, since they all contribute to the coefficients of the scalar D-functions occurring in 
the one-loop Zbb amplitude. For instance, it has been particularly useful in the case of box 
diagrams like the one shown in Fig. [3j since this diagram and related ones contain up to 
DA-PV functions that cannot be reduced even at the level of the amplitude squared. Since 
they involve up to four powers of inverse GDs, they are particularly subject to numerical 
instabilities and it is important to have their analytic expressions as compact as possible. 

After the tensor integral reduction is performed, the fundamental building blocks are 
one-loop scalar integrals with up to five denominators. They may be finite or contain both 
ultraviolet (UV) and infrared (IR) divergences. The finite scalar integrals are evaluated using 
the method described in Ref. 52j and cross-checked with the numerical package FF 44]. 
The UV and IR singular scalar integrals are calculated analytically by using dimensional 
regularization in d = 4 — 2e dimensions. The most difficult integrals arise from IR divergent 
pentagon diagrams with several external and internal massive particles. We calculate them 
as linear combinations of box integrals using the method of Refs. 55|, |56( and of Ref. 52]. 
Details of the box scalar integrals (see also Ref. Q) and the pentagon reduction, as well as 
the set of IR-divergent three and two-point functions used in this calculation, are given in 
Ref. [50] . 



The UV singularities of the virtual cross section are removed by introducing a suitable 
set of counterterms (see Refs. 48, S] for details), while the residual renormalization 
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scale dependence is checked from first principles using renormalization group arguments as 
in Eq. (4) of Ref. 49j. Note that we use the on-shell subtraction scheme when fixing the 
wave function renormalization constant of the external bottom quark field (SZ^) and the 
mass renormalization constant (Srrif,). The IR singularities of the virtual cross section are 
canceled by analogous singularities in the 0(aa^) real cross section. 

In our calculation we treat 75 according to the naive dimensional regularization approach, 
i.e. we enforce the fact that 75 anticommutes with all other 7 matrices in d = 4 — 2e 
dimensions. This is known to give rise to inconsistencies when, at the same time, the d- 
dimensional trace of four 7 matrices and one 75 is forced to be non-zero (as in d = 4, 
where Tr(7 At 7 I/ 7 p 7°"7 5 ) = Aie^ vpu ) [581 ] . In our calculation, both UV and IR divergences are 
handled in such a way that we never have to enforce simultaneously these two properties 
of the Dirac algebra in d dimensions. For instance, the UV divergences are extracted and 
canceled at the amplitude level, after which the d — > 4 limit is taken and the renormalized 
amplitude is squared using d — 4. Thus, all fermion traces appearing at this point are 
computed in four dimensions and therefore have no ambiguities. 

We note that the tree-level amplitude Aq in Eq. ([7]) has generically to be considered as a d- 



dimensional tree-level amplitude. This matters when the J- 
IR divergent. Actually, as it has been shown in Refs. 48 



n amplitudes in Eq. ([7]) are UV or 



49l | . both UV and IR divergences 



are always proportional to the tree level amplitudes and they can be formally canceled 
without having to explicitly specify the dimensionality of the tree level amplitude itself. 
After UV and IR singularities have been canceled, the remaining phase space integration is 
computed in d — 4 dimensions using standard MC techniques. 



C. The real cross section (rj| al 

The NLO QCD real cross section <5jJ al in Eq. ([3]) corresponds to the 0(a s ) corrections to 
ij — > Zbb due to the emission of an additional real parton, i.e. to the process ij — > Zbb + g, 
and the tree-level process q(q)g — > Zbb+q(q). <r[J contains IR singularities which cancel the 
analo gou s singularities present in the 0(a s ) virtual corrections and in the NLO PDFs (see 



Refs. 



48 



49 



50| for details). These singularities can be either soft, when the emitted extra 
parton is a gluon and its energy becomes very small, or collinear, when the final-state parton 
is emitted collinear to one of the partons in the initial state. There is no collinear singularity 



10 



arising from the radiation off the final-state bottom quarks, since they are considered to be 
massive. 

We have calculated the cross sections for the processes 

%i) + jfe) -> KPb) + KPb) + z (Pz) + g(k) 

and 

(Q, Q)(Qi) + 9(<h) ~> KPb) + KPl) + Z(p z ) + (q, q)(k) , 
with qi + q 2 = Pb + Pi + Pz + k, using the two-cutoff Phase Space Slicing (PSS) method. 
This implementation of the PSS method was originally developed to study QCD corrections 
to dihadron production |59] and has since then been applied to a variety of processes (a nice 



review can be found in Ref. 



60|). We fo 



low closely the application of the PSS method to 
to which we refer for more extensive references 



Htt production as presented in Refs. [4£ 
and full details. Although we are considering Zbb production, the kinematics are equivalent, 
and the color structure and IR behavior are the same, so necessarily their soft and collinear 
kernels are the same. In the following we briefly summarize our implementation of the 
two-cutoff PSS method. 

Using the PSS method, the IR singularities can be conveniently isolated by slicing the 
phase space of the final-state particles into different regions defined by suitable cutoffs. To 
isolate the soft and collinear singularities we impose soft (S s ) and collinear (S c ) cutoffs on the 
phase space of the emitted parton as follows. By introducing an arbitrary small soft cutoff 
8 S , we separate the overall integration of the qq, gg — > bbZ + g phase space into two regions 
according to whether the energy of the final state gluon (k° = E g ) is soft, i.e. E g < 6 s s/s/2, 
or hard, i.e. E g > 5 s y/s/2. In order to isolate the collinear singularities we further divide 
the hard region of the qq,gg — > bbZ + g phase space into a hard/ collinear and a hard/non- 
collinear region, by introducing a second small collinear cutoff 5 C . The hard/non- collinear 
region is defined by the condition that both 

J^ 1 ■= > 5 C and J? 2 = > 5 C (8) 

are true. We apply the same collinear cutoff to the tree- level process q(q)g — > Zbb + q(q). 

rpi i j iv , r 4-l_ 1 i.- * hard/non— coll c j i 

lhe hard non-collmear parts of the real cross sections, o~ q q gg qg , are finite and can be 
computed numerically. The partonic real cross sections can then be written as follows: 

^real _ '•soft . a hard/coll . a hard/non— coll /q\ 

a qq,gg,qg ~ a qq,gg + °~qq,gg,qg + °~qq,gg,qg > W ) 

11 



where crfj^ and v^gg^ 1 is obtained by integrating analytically over the soft and collinear 
regions of the phase space of the emitted parton, respectively, and contains all the IR 
divergences of cr qq ^ ggqg . The dependence on these arbitrary cutoffs, d~ s , 5 C , is not physical, 
and cancels at the level of the real cross section, i.e. in <jjf al . This cancellation constitutes 
an important check of the calculation. 

We conclude this section by showing explicitly that the total hadronic cross section at 
NLO QCD does not depend on the arbitrary cutoffs introduced by the PSS method, i.e. on 
5 S and 5 C . The cancellation of the PSS cutoff dependence is realized in o"[| al by matching 
contributions that are calculated either analytically (cr^^ and frhard/coii^ ^ e iR-unsafe re- 
gion below the cutoffs, or numerically, in the IR-safe region above the cutoffs (jj h ^ rd l non coll y 
While the analytical calculation in the IR-unsafe region reproduces the form of the cross 
section in the soft or collinear limits and is therefore only accurate for small values of the 
cutoffs, the numerical integration in the IR-safe region becomes unstable for very small val- 
ues of the cutoffs. Therefore, obtaining a convincing cutoff independence involves a delicate 
balance between the previous antagonistic requirements and ultimately dictates the choice 
of values that are neither too large nor too small for the cutoffs. In Figs. H] and El using 
the setup described in Section HTT1 we demonstrate the independence of a NLO (pp — > Zbb) on 
5 S and 5 C separately, by varying only one of the two cutoffs over an extended range, while 
the other is kept fixed. In Fig. HI d~ s is varied between 10~ 5 and 1CT 2 with 5 C = 10~ 5 , while 
in Fig. El o~c is varied between 10~ 7 and 10~ 4 with 5 S = 1CT 3 . In both plots, we show in the 
upper window the overall cutoff dependence cancellation between the hadronic cross sections 

E, soft , hard/coll\ -i v—^ hard /non— coll . rpa i . , -,. u 1 1 — 1 

iji^ij + <7ij ) and l^ij a ij m l^ij a ij ' including all channels, gg, qq and 

qg. Note that we also take into account contributions from the LO and the virtual cross sec- 
tions which are cutoff independent. In the lower window of the same plots we show the full 
<7 NLO , including all channels, on a scale that magnifies the details of the cutoff-dependence 
cancellation. The statistical errors from the MC phase space integration are also shown. 
Both Figs. H] and El show a clear plateau over a wide range of S s and 5 C and the NLO cross 
section is proven to be cutoff independent. The numerical results presented in Section II II I 
have been obtained by using the two-cutoff PSS method with 5 S = 10~ 3 and <5 c = lCr 5 . 
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FIG. 4: Dependence of a NLO (pp — > Zbb) on the soft cutoff 5 S of the two-cutoff PSS method for 
fi = 2m b + M z , and<5 c = 10- 5 . The upper plot shows the cancellation of the ^-dependence between 
(jsoft _j_ ^hard/coll anc j ^hard/non-coii _ rj^g j ower pj t shows, on an enlarged scale, the dependence 

of the full a NLO 
integration. 
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III. NUMERICAL RESULTS 

The results for Zbb observables presented in this paper are obtained for the Tevatron 
pp collider at sh = 1-96 TeV. If not stated otherwise, we assume a non-zero bottom-quark 
mass, fixed at mj = 4.62 GeV. The mass of the top quark, entering in the virtual corrections, 



is set to m t = 170.9 GeV. The Z-boson mass is taken to be Mz = 91.1876 GeV [6lJ and 
the IV-boson mass is calculated from Mw = Mz cos9 w with sin 2 ^ = 0.223. We work in 
the electroweak G M input scheme and replace the fine structure constant a(0) = e 2 /(47r) 
by ot,(Gfj) = ^G^M^ sin 2 6 W with the Fermi constant G M = 1.16639 ■ 10- 5 GeV" 2 . The 



LO results use the one-loop evolution of a s and the CTEQ6L1 set of PDFs 62], with 
a^°(Mz) = 0.130, while the NLO results use the two- loop evolution of a s and the CTEQ6M 
set of PDFs, with a^ LO (Mz) = 0.118. In the calculation of the parton luminosity we assume 



13 



40 
30 
20 

I 10 

O 

e" 
-10 
-20 



3.44 
3.4 
3.36 
3.32 



H 1 1 1 — I I I I | 1 1 1 1 — I I I I | 1 1 1 1 — I I I- 

Total (all channels) NLO Inclusive 

-- Soft+Coll+Virt+Tree 
Hard non-coil 



5 = 10 

S 

^. = ^. = M z + 2m b 



2->4 



cuts: p ( > 15 GeV 
N<2 
R = 0.7 

H 1 I I I I I 1 1 



H 1 I I I I I I 



2->3 

H 1 I I I I I I 



_i i i i i i i i i i i i i i i i i i i i i i i i i i 



le-07 



le-06 



le-05 



0.0001 



FIG. 5: Dependence of a Nh °(pp — ► Zbb) on the collinear cutoff 5 C of the two-cutoff PSS method, for 
H = 2m b + M z , and <5 S = 1CT 3 . The upper plot shows the cancellation of the 5 s -dependence between 
a soft _|_ a hard/coii ^ anc ? a hard/non-coii _ rpj^g i ower ^\ ^ s h OWS) on an enlarged scale, the dependence 

of the full <t nlo = (Tgg ° + o"^ + a qg° on with the corresponding statistical errors of the MC 



integration. 



five light flavors in the initial state. Including the 6-quark PDF has a negligible effect 



< 0.1%) on the Zbb cross section and is inc 



implement the jet algorithm 63 



64 



65 



uded to consistently compare with MCFM. We 



66j with a pseudo-cone size R = 0.7 and we 



recombine the parton momenta within a jet using the so called covariant i?-scheme [64j. We 
checked that our implementation of the kx jet algorithm coincides with the one in MCFM. 
We require all events to have a bb jet pair in the final state, with a transverse momentum 
larger than 15 GeV (p^ b > 15 GeV) and a pseudorapidity that satisfies \r] b '^\ < 2. We impose 
the same pt and \r]\ cuts also on the extra jet that may arise due to hard non-collinear real 
emission of a parton, i.e. in the processes Zbb + g or Zbb + q(q). This hard non-collinear 
extra parton is treated either inclusively or exclusively. In the inclusive case we include 
both two- and three-jet events, while in the exclusive case we require exactly two jets in the 
event. Two-jet events consist of a bottom-quark jet pair that may also include a final-state 
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light parton (gluon or quark) due to the applied recombination procedure. Results in the 



30] 



massless bottom-quark approximation have been obtained using the MCFM code 

In Table |T] we present results for the total LO and NLO QCD pp — > Zbb cross sections, 
obtained with the scale \x r = fif = Mz+2mb, for both our fully massive calculation and in the 
massless approximation. As can be seen, the NLO QCD corrections increase considerably 



TABLE I: LO and NLO total Zbb cross sections at the Tevatron for massive and massless bottom 
quarks, using \i r = fj,f = Mz + 2m&. The numbers in square brackets are the ratios of the NLO 
and LO cross sections, the so called K-factors. Statistical errors of the MC integration amount to 
about 0.1%. 



Cross 


Section 


m b (pb) [ratio] m& 


= (pb) [ratio] 






2.21[-] 


2.37[-] 


<x NLO 


inclusive 


3.40[1.54] 


3.64[1.54] 


a NLO 


exclusive 


2.80[1.27] 


3.01[1.27] 



the total cross section, with NLO vs. LO ratios (fT-factors) that, in both the massive 
and massless bottom-quark case, amount to K — 1.54 and K = 1.27 for the inclusive 
and exclusive case, respectively. In the following we will study the impact of the NLO 
QCD corrections on Zbb observables in more detail. Specifically we will show examples of 
kinematic distributions where a global rescaling (or i^-factor) does not properly describe the 
effect of these corrections. In Figs. [6] and [7] we illustrate the renormalization and factorization 
scale dependence of the LO and NLO QCD total cross sections, both in the inclusive and 
exclusive case. Figure [6] shows the overall scale dependence of both LO, NLO inclusive and 
NLO exclusive total cross sections, when both fi r and /// are varied independently between 
Ho/2 and 4// (with /x = m b + M z /2), including full bottom-quark mass effects. We notice 
that the NLO QCD cross sections have a reduced scale dependence over the range of scales 
shown, and the exclusive NLO QCD cross section is more stable than the inclusive one. 
This effect is mainly driven by the tree level subprocess q(q)g — > Zbb + q(q) contributing 



to the real corrections. This is illustrated by the right hand side (RHS) plots of Figs. 7(a) 



and 7(b) , where we show separately the /^-dependence of the total cross sections to the qq, 
qg + qg and gg initiated processes, for /i r = [if, both for the inclusive and for the exclusive 
case. It is clear that the low scale behavior of the inclusive cross section is considerably 
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FIG. 6: Dependence of the LO (black solid band), NLO inclusive (blue dashed band), and NLO 
exclusive (red dotted band) Zbb total cross sections on the renormalization/factorization scales, 
including full bottom-quark mass effects. The bands are obtained by independently varying both 
fjL r and /j,f between uq/2 and 4//q (with /jlq = ra^ + Mz/2). 



affected by the qg + qg contribution, which show a monotonic dependence on [i (i.e. with no 



plateau) characteristic of tree level processes. In the left hand side (LHS) plots of Figs. 7(a 



and 7(b) we also compare the scale dependence of our results to the scale dependence of the 
corresponding results obtained with = (using MCFM), both at LO and at NLO. Using 
a non-zero value of is expected to have a small impact on the scale dependence of the 
results 2 , since the only modification to the renormalization scale dependence originates from 

n 

the bottom-quark mass and field renormalization, as discussed in Section IIB of Ref. |27| . 
where we compare the minimal and on-shell subtraction schemes. Indeed, as can be seen 



in Figs. 7(a) , 7(b) the scale dependence of the LO and NLO curves is very similar for both 
the case of a massive and massless bottom quark. While the LO cross section still has 
a 45% uncertainty due to scale dependence, this uncertainty is reduced at NLO to about 
20% for the inclusive and to about 11% for the exclusive cross sections. The uncertainties 
have been estimated as the positive/negative deviation with respect to the mid-point of the 
bands plotted in Fig. (6J where each band range is denned by the minimum and maximum 



2 Note that we always use m& = 4.62 GeV in the determination of the scales in terms of uq = nib + Mz/2 
even in the results obtained with mb = 0. 
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(a,) Inclusive case 




(b) Exclusive case 



FIG. 7: Dependence of the LO and NLO inclusive and exclusive Zbb total cross section on the 
renormalization/factorization scale, when \i r = fj,f = /x. The LHS plots compare both LO and 
NLO total cross sections for the case in which the bottom quark is treated as massless (MCFM) 
or massive (our calculation). The RHS plots show separately, for the massive case only, the scale 
dependence of the qq, gg and qg + qg contributions, as well as their sum. 



17 



0.06 - h> = M z /2 + m b 



0.03 



• Inclusive case 
A Esclusive case 



NLO NLO , LO , LO „ 





< 



-0.03 



-0.06 



0.5 



cuts: p ( > 15 GeV 

mi < 2 

R = 0.7 

_J I 



FIG. 8: Dependence on the renormalization/factorization scale of the rescaled difference between 
our NLO calculation (with mj 7^ 0) of the Zbb total cross section and MCFM (with m b = 0) for 
the inclusive and exclusive cases (with /x r = ///). The error bars indicate the statistical uncertainty 
of the MC integration. 



value in the band. We notice incidentally that the difference in the total cross section due 
to finite bottom-quark mass effects is less significant than the theoretical uncertainty due to 
the residual scale dependence in the inclusive case, but is comparable in size in the exclusive 
case. Indeed, the finite bottom-quark mass effects amount to a reduction of the total cross 
sections by about 7% compared to the massless case at both LO and NLO QCD. 

In Fig. [8], we show the rescaled difference between the NLO total cross sections obtained 
from our calculation (with m b 7^ 0) and with MCFM (with m;, = 0) denned as follows: 

Aa = a NLO (m b ? 0) - a^°(m b = 0) ± °> . (10) 

a LU (m b = 0) 

As can be seen, within the statistical errors of the MC integration, the finite bottom-quark 
mass effects on the total cross sections at NLO are well described by the corresponding 
effects at LO. 

Finally, in Figs.[H]to[TT]we study the distribution da/dm b i, where m bb is the invariant mass 
of the bb jet pair. The impact of NLO QCD corrections on this distribution is illustrated in 



Figs. 9(a) and 9(b) for the inclusive and exclusive case, respectively. We see that the NLO 



QCD corrections affect the differential cross section quite substantially. In each figure the 
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FIG. 9: The distribution da/dm bb at LO and NLO QCD. The RHS plots show the ratio of the LO 
and NLO distributions. 

RHS plot gives the ratio of the NLO and LO distributions. We stress the fact that the NLO 
m bb distributions cannot be obtained from the LO ones by just rescaling, which is clear from 
the RHS plots of Fig. E 



Figs. 10(a) and 10(b) compare the NLO do jdm bb distributions obtained from the massive 
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FIG. 10: The inclusive and exclusive NLO QCD distributions dcr/dm bb - derived from our calculation 
(with nif) ^ 0) and from MCFM (with nib = 0). The RHS plots show the ratio of the two 
distributions, da(mh ^ 0)/ da (mi, = 0). 
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and massless bottom-quark calculations. The results with m b = have been obtained using 
MCFM. As expected, most of the difference between the massless and massive bottom- 
quark cross sections is coming from the region of low m bb invariant mass, both for the 
inclusive and exclusive case, where the cross sections for mj 7^ are consistently below 
the ones with m;, = 0. This is emphasized in the RHS plots, where we show the ratio of 
the two distributions, da{m b 7^ 0)/da(m b = 0). For completeness, we also show in Fig. [11] 
the comparison between massive (m b 7^ 0) and massless (nib = 0) calculations at LO in 
QCD. The LO m b i distribution for massive bottom-quarks has been obtained both from our 
calculation and from MCFM, which implements the m;, 7^ option at tree level, and both 
results agree perfectly. In general, mass effects are similar at LO and NLO. To illustrate 
this in more detail we show in Fig. [T5] the rescaled difference between the m h i distributions 
obtained with our NLO calculation (with m;, 7^ 0) and with MCFM (with m& = 0) defined 
as follows: 

da rfcx NLO , , , rfa NLO , ,da hO (m b ^0) , , 

= (m b ^ (m b = , TM T . 11 

dm bb dm bb dm bb d<j LU (m b = 0) 

We notice that, in the Zbb case, finite bottom-quark mass effects are relevant up to values 

of the m bb invariant mass around 50 GeV. Although not included in the present analysis, 

our calculation is still valid when both b quarks are in the forward direction. In this region, 

collinear singularities can arise, which are regularized by the finite 6-quark mass. The 

resummation of the corresponding large logarithms is then appropriate and is left to future 

improvements. 



IV. CONCLUSIONS 

We have calculated the NLO QCD corrections to hadronic Zbb production including 
full bottom-quark mass effects. We have presented numerical results for the total cross 
section and the invariant mass distribution of the bottom-quark pair at the Tevatron for 
both massless (with MCFM) and massive bottom quarks. We apply the kx jet algorithm, 
require two 6-tagged jets and impose kinematical cuts that are inspired by the D0 and CDF 
searches for the SM Higgs boson in ZH production. The NLO QCD Zbb cross section shows 
a considerably reduced renormalization and factorization scale dependence, i. e. about 20% 
for the inclusive and about 11% for the exclusive cross sections as opposed to 45% scale 
uncertainty of the LO cross section. The bottom-quark mass effects amount to about 7% 
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FIG. 11: The LO distribution da/dm b i derived from our calculation (with m& 7^ 0) and from MCFM 
(with m& = 0). The RHS plot shows the ratio of the two distributions, da (nib 7^ 0)/ do (nib = 0). 
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FIG. 12: The m b i distribution of the rescaled difference between our NLO calculation (with nib 7^ 0) 
and MCFM (with nib = 0) f° r the inclusive (upper plot) and exclusive case (lower plot). 
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of the total NLO QCD cross section and can impact the shape of the m h i distributions, in 
particular in regions of low m b i. This is relevant to SM Higgs searches in ZH associated 
production with H — > bb and to searches for MSSM Higgs bosons in Hbb production with 
H — > fi + fi~,T + r~. We also plan to apply the formalism developed in this paper to the 



calculation of both Ztt 



671 . l68j | and ■ytt production at NLO in QCD. Both processes are 

m, |7q|, while Ztt also 



of interest to the study of electroweak properties of the top quark 
constitutes a relevant background to new physics searches. 
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